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Abstract 


We give a determination of the phenomenological value of the Wilson 
(or gradient) flow scales Iq and wq for 2 + 1 flavours of dynamical quarks. 
The simulations are performed keeping the average quark mass constant, 
which allows the approach to the physical point to be made in a controlled 
manner. 0{a) improved clover fermions are used and together with four 
lattice spacings this allows the continuum extrapolation to be taken. 


1 Introduction 


Numerical lattice QCD simulations naturally determine dimensionless quantities 
such as mass ratios and matrix element ratios, however determining a physical 
value requires the introduction of a scale, usually taken from experiment. A 
hadron mass, such as the proton mass, or decay constant, such as the pion decay 
constant, are often used for this purpose. We discuss here setting the scale 
using flavour-singlet quantities, which in conjunction with simulations keeping 
the average quark mass constant allow SU{3) flavour breaking expansions to be 
used. This is illustrated here using 2 + 1 clover fermions, and a determination 
of the Wilson flow scales to and Wq is given. These are ‘secondary’ scales and 
are not experimentally accessible and thus they have to be matched to physical 
quantities. These flow scales are cheap to compute from lattice simulations (for 
example they do not require a knowledge of quark propagators) and accurate (for 
example they do not require a determination of the potential which requires the 
limit of a large distance). So once the phenomenological value of the flow scales 
is known the determination of physical values becomes more tractable. 

Flow and flow variables were introduced by Liischer, [1]. We follow him 
here, [2], in particular in our brief discussion of the to scale. Flow represents 
a smoothing of the gauge helds. We denote the flow time by t, and the link 
variables at this time by U^{x,t) = exp(iT“6*“(a:, t)) which evolve according to 


dUf,{x,t) 

dt 


5ei{x,t) 


u^{x,t) , 


with f/^(a;,0) 


U,{x) , 


( 1 ) 


with S'flow[f^] being the flow action, which does not have to be the same as the 
action used to generate the gauge variable, (x is just the normal 4-dimensional 
Euclidean space-time.) Setting 


F(t)=f{E(t)), where £(4) =++«), (2) 

then we define the to scale by 


^WIt=to(c)=C. 


( 3 ) 
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Alternatively, [3] define the wq scale as 



= c, 


(4) 


where in both dehnitions c is a constant, conventionally taken as c = 0.3. We 
require a value of c such that a \/8io L, L being the lattice size and this 
value was found to be a suitable choice, [2]. Alternative suggestions have been 
made, see e.g. mu- 



In this article we shall give a determination of J and There have been 


[m 121 [ig. We have also published preliminary results, [IT] . 

The plan of this article is as follows. In section |2] we describe our method of 
approaching the physical quark mass starting from a point on the SU (3) flavour 
symmetric line, [13 US]- We also discuss the general property of singlet quanti¬ 
ties, that they have a stationary point about this SU{3) flavour symmetric line. 
Section [3] gives examples of singlet quantities both hadronic and gluonic (i.e. in 
this case to and Wq) and also discusses their SU (3) flavour breaking expansions. 
Section m first gives our lattice conventions, ensembles used and numerical values 
of the singlet quantities. This is followed by section [5] in which the i/fo and wq 
scales are determined for several lattice spacings. In the next section, section Owe 
take the continuum result to give the final result. Finally in section [3 we compare 
our result with other results for rif = 2 + 1 flavours and give our conclusions. 

2 Extrapolating flavour singlet quantities 

We consider extrapolations to the physical point from a point on the SU{3) 
flavour symmetric line keeping the average quark mass fixed, [iac6], 


m = ^{mu + rrid -|- mg) = const.. 


(5) 


This means that as the pion mass tends downwards to its physical value, the 
kaon mass increases upwards to its physical value. (In particular the kaon mass 
is never larger than its physical value.) Possible scenarios are sketched in Fig. [1] 
for a path from a point on the SU (3) flavour symmetric line to the physical point 
for the case of mass degenerate u and d quarks mi, see eq. flT^) . (For non-mass 
degenerate u and d quarks we would have instead a plane.) Shown in Fig. [1] 
are the bare and renormalised quark masses for the case discussed here of clover 
fermions. (Because the singlet and non-singlet pieces renormalise differently the 
renormalised quark axes are not orthogonal to each other, as further discussed in 
|16] . For chiral fermions this would not be the case.) In the left hand panel one 
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Figure 1: Possible scenarios for the path in the (m„ = =) mi - mg plane. In the 

LH panel m = const., while in the RH panel we hold other singlet quantities constant. 


common trajectory is sufficient, while in the right hand panel it depends on the 
singlet quantity used, [T6] . 

With the condition of eq. ([5]) flavour singlet quantities turn out to have a sta¬ 
tionary point in the quark mass starting from a given point on the SU{3) flavour 
symmetric line. As we shall see this potentially allows simpler extrapolations to 
the physical point. This property may be shown by considering small changes 
about a given point on the SU{3) flavour symmetric line. Let Xs{mu,md,ms) 
be a flavour singlet object i.e. Xs is invariant under the quark permutation sym¬ 
metry between u, d and s. Then Taylor expanding Xs about a point mo on the 
SU{3) flavour symmetric line m„ = m^ = rris = mo = rn, 

niq = m + 6mq , (6) 


gives 


Xs{m + Sniui m -|- 6md, m + 6ms) 


= Xs{m, m, m) -f 


dXs 


dm,, 


6mu + 


dXs 


dm„ 


6md + 


dXs 

dmg 


But on the symmetric line we have 


dXs 

dXs 

dXs 

dmu 

1 

o 

0 dms 


and on our chosen trajectory, eq. (j5|) 

6mu + 6md + 6ms = 0 , 


(7) 

6ms + 0{{6mq)'^). 

( 8 ) 

(9) 
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which together imply that 

Xsijn + m + 6md, m + 5ms) = Xs(m, m, m) + 0{{5mq)‘^). ( 10 ) 

In other words, the effect at first order of changing the strange quark mass is 
cancelled by the change in the light quark mass, so we know that Xs must have 
a stationary point on the SU{3) flavour symmetric line. 

Thus provided the quadratic terms in eq. oa are small it will not matter 
which singlet quantity we use, all are equivalent to keeping m = const., the 
scenario of the left hand panel of Fig. [TJ 

3 Defining singlet quantities 

We shall consider here hadronic singlet quantities, such as independent aver¬ 
ages with respect to quark permutations of the pseudoscalar, vector and nucleon 
octets, [ini [16], and also gluonic quantities derived from the Wilson flow, pis]. 
In particular we shall consider 

Xl = |(M|.,+M|o + M2++M2_+M|._+M|o) = i(2Mi + M2) 

+ M2+ + M2_ + M|,o + M|._) = + M^) 

Xjf = \{Ml + MI + + M|_ + M|o + M|_) = \{M% + M| + Ml) 

K = i/to 

^^0 = IK- (11) 

(The charge conjugate mesons have the same masses, at least for pure QCD. For 
example M^+ = M^--, Mk+ = and M^o = T/^o, obviously no such similar 
result holds for the baryons.) The second expressions are the masses for mass 
degenerate u and d quarks, 

mu = md = mi. ( 12 ) 

Except for X^, which is naturally an average over quadratic masses, it does not 
matter whether we consider linear or quadratic averages - quadratic averages 
were found to give slightly better fits for heavy partially quenched masses (up 
to the charm quark); in the small quark mass range considered here this is less 
important. Note also that for the Wilson flow singlets we consider the inverse of 
the flow variable as then all singlet quantities have the same dimensions. 

Other possibilities, as discussed in |T6] include the further nucleon octet singlet 

Xl = l(Mlo + M|o) = i(Ml + Mi), (13) 

and baryon decuplet singlets 

Xi = \{Ml^^+Ml.+Ml.) = \{2Ml + Ml) 

XL = |(Mi++Mio + M2.+ +M2._+M|*o + M|._) = 2(Mi + M| + M|) 

X|* = M|.o = ML . (14) 
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Singlet qnantity 

GeV 

Xexp 

0.4126 

x^^p 

0.8562 

\^exp 

1.1610 

j^exp 

1.1548 

j^exp 

1.3944(12) 

j^exp 

1.3837(1) 

-V^exp 

1.3888(6) 


Table 1: Experimental values for various singlet quantities, averaging over (masses)^ 
(to 4 decimal places). 


Further possibilities can be constructed from ‘fictitious’ particles, such as a ‘nu¬ 
cleon’, Ng, with three mass degenerate quarks at the strange quark mass. At the 
SU{3) flavour symmetric points, all the baryon octet hadrons are mass degenerate 
(at least from a QCD perspective), so we expect that away from this point there 
is no (or very little) difference between Xff and i.e. = X\. The same ar¬ 
gument holds for the various baryon decuplet possibilities, i.e. Xa = = X^*. 

How far does this extend? Let us consider the experimental (or phenomeno¬ 
logical) singlet hadron mass result^ as given in Table [^. It is seen that even 
after the extrapolation from the SU{3) flavour line to the experimental point, 
then X“^ ^ X^'’ and X“^ ^ XS’’ X^^ the worst discrepancy (between XS'’ 
and X|;?) is only a fraction of a percent. This indicates that quite likely the X 5 
are constant over a large interval. 

However as the baryon decuplet possibilities are numerically substantially 
noisier, we shall only use the baryon octet hadrons here. 

An equivalent statement (for singlet quantities built from hadron masses) is 
found by considering the SU{3) flavour breaking expansion. As discussed in 
usms] we have for the octet mesons, n~^{ud), vr (du), X+(ms), K (su), K^{ds) 
and {sd) not lying at the centre of the octet 

M‘^{ab) = Mq 4- aTr{Sma + drrih) + 0((5mg)^), (15) 

where a and b are u, dor s quarks. Similar results hold for the vector mesons. For 
the p{uud), n{ddu), S+(mms), T,~{dds), E^{ssu), E~{ssd) baryons on the outer 

^Only when necessary and for clarity do we distinguish between experimental and lattice 
masses - and Xlf* respectively. 

^As we are not considering mass differences, then the effect of electromagnetic effects is small, 
and so we can disregard them here. For example for the lightest particles - the pseudoscalar 
octet, the value given in Table [T] for Ar,r is to be compared with the value upon using Dashen’s 
theorem which gives 0.4116GeV (see e.g. [E]). This is a < 0.2% difference. For the baryons 
for Xn it is < 0.1% 


6 









ring of the baryon octet0 we have 

M‘^{aab) = Mq^ + Ai{2Sma + Snib) + A 2 {Smb - 5ma) + 0{{5mqY). (16) 


All the expansion coefficients are functions of m. It is easy to check that this 
means that X| = + 0((5m/)^), S = ti, p and N in agreement with eq. (1T0|1 . 

‘Fan’ plots from the symmetric point down to the physical point are well 
described by the linear behaviour of eqs. (|T5il and (|T6|l as shown in [151 CHI [IHl EH] 
which further supports the earlier statement that Xs is constant over a large 
quark mass range. 

Although in our extrapolations, we shall not be using chiral perturbation 
theory, yPT, it is natural to ask about its relationship to the SU{3) flavour 
breaking expansion. This also allows a check on eq. (TO. assuming we are in 
a region where yPT is valid. This was investigated in [T5l 116] for hadron mass 
singlets and we now extend the argument to to (and Wq), using the result of [21]. 
Using the notation of this paper and for mass degenerate u and d quark masses 
we hnd 


to = T{x) 


1 + 


(4vr/o) 


,{l^2 + \k'l){Xs-Xi? + 


where 


T{x) = to,ch 


1 + 


3fci 


■x + 


8 kb 


—2 1 X I 


9k' 


X 


(17) 


( 18 ) 


(4!r/o)‘^ (4!r/o)* (lir/o)"* 

being the value of to on the symmetric line (to,ch is the value in the chiral limit), ki 
are constants, fo the pion decay constant again in the chiral limit and xi = Bomi, 
Xs = Boms, X = |(2xz + Xs)- 

As expected, there is no linear term, and the hrst term we see is quadratic in 
the SU{3) breaking. Further details are given in Appendix lAl 


4 Lattice matters 

4.1 General 

We consider 2 + 1 non-perturbatively 0(a) improved clover fermions, as described 
in ca. The relation between the bare quark masses in lattice units and the 
lattice mass parameters Kg is given by [laiis] 

mg = ^ -^ I with q G {/, s, 0} , (19) 

2 y Kg KOc ) 

•^For non-degenerate u and d quark mass, the Lambda and Sigma particles mix. This mixing 
is however very small. This was investigated in |20] . where using the notation there, it was 
shown that ^(M^o + M^o) = Pa^ = +toAr + 0((5mqY), i.e. again we have no linear term in 
the quark mass. Note also that when the u and d quark masses are degenerate, then no mixing 
occurs. 
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/S 

V 

Ko 

Hi 

Hs 

M^L 

L [fm] 

5.80 

48^ X 96 

0.122760 

0.122760 

0.122760 

6.95 

2.82 

5.80 

48^ X 96 

0.122810 

0.122810 

0.122810 

6.11 

2.82 

5.80 

48^ X 96 

0.122810 

0.122880 

0.122670 

5.11 

2.82 

5.80 

48^ X 96 

0.122810 

0.122940 

0.122551 

4.01 

2.82 

5.80 

48^ X 96 

0.122870 

0.122870 

0.122870 

4.96 

2.82 

Table 2: Parameters for (3 = 5.80. 

Each block has the same kq, i.e. constant m 

/S 

V 

Ko 

Hi 

Kg 

M^L 

L [fm] 

5.65 

32^ X 64 

0.121975 

0.121975 

0.121975 

4.99 

2.19 

5.65 

32^ X 64 

0.122005 

0.122005 

0.122005 

4.67 

2.19 

5.65 

323 X g4 

0.122005 

0.122078 

0.121859 

4.00 

2.19 

5.65 

323 X g4 

0.122005 

0.122130 

0.121756 

344 

2.19 

5.65 

323 X g4 

0.122030 

0.122030 

0.122030 

4.32 

2.19 

5.65 

323 X g4 

0.122050 

0.122050 

0.122050 

4.03 

2.19 


Table 3: Parameters for (3 = 5.65. The entries in italics have < 4. 


where in simulations we have mass degenerate u and d quarks, i.e. = ma = mi 
and the s quark has mass m^. Along the SU{3) flavour mass degenerate line, the 
common quark mass is denoted by mo (or equivalently by kq) and where vanishing 
of the quark mass along this line determines kqc- Along the m = mo = constant 
line gives from eqs. ([6]) and flT^ the SU{3) flavour breaking mass parameter as 

f ^ 

2\Kg KoJ 

We see that Koc has dropped out of this equation, so we do not need its explicit 
value here. Along this trajectory the choice of quark masses is restricted and we 
have 

1 , , 
~ _3_g ’ 

Ko Kl 

SO once we have decided on a kq, then a given ki determines Kg. 

We consider four beta values [3 = 5.8, 5.65, 5.50, 5.40 (where (3 = lO/f^g with 
our conventions). In Tables [21 [3], H] and [5] we give parameters of the runs. The 
entries in italics have < 4. 

For example we consider f3 = 5.50 on 32^ x 64 lattices with degenerate quark 
masses of kq = 0.120900, 0.120920, 0.120950 and 0.120990 which, as we will 
see, encompasses the initial SU{3) flavour symmetric point on the constant m 
trajectory to the physical point. As can be seen for some of the kq values we have 


( 20 ) 














/S 

V 

Ko 


Kg 

M^L 

L [fm] 

5.50 

32^^ X 64 

0.120900 

0.120900 

0.120900 

5.59 

2.37 

5.50 

32^ X 64 

0.120900 

0.121040 

0.120620 

4.32 

2.37 

5.50 

323 X g4 

0.120900 

0.121095 

0.120512 

3.72 

2.37 

5.50 

323 X g4 

0.120900 

0.121145 

0.120413 

3.10 

2.37 

5.50 

483 X 96 

0.120900 

0.121166 

0.120371 

4.10 

3.55 

5.50 

323 X g4 

0.120920 

0.120920 

0.120920 

5.27 

2.37 

5.50 

323 X g4 

0.120920 

0.121050 

0.120661 

4.10 

2.37 

5.50 

323 ^ g4 

0.120950 

0.120950 

0.120950 

4.83 

2.37 

5.50 

323 X g4 

0.120950 

0.121040 

0.120770 

3.97 

2.37 

5.50 

323 X g4 

0.120950 

0.121099 

0.120653 

3.24 

2.37 

5.50 

323 X g4 

0.120990 

0.120990 

0.120990 

4.11 

2.37 


Table 4: Parameters for (3 = 5.50. 



V 

Ko 


Kg 

M^L 

L[fm] 

5.40 

243 X 48 

0.119860 

0.119860 

0.119860 

4.98 

1.96 

5.40 

243 X 48 

0.119895 

0.119895 

0.119895 

4.54 

1.96 

5.40 

243 X 48 

0.119930 

0.119930 

0.119930 

4.11 

1.96 

5.40 

243 X 48 

0.119930 

0.120048 

0.119695 

3.35 

1.96 

5.40 

243 X 48 

0.120000 

0.120000 

0.120000 

3.25 

1.96 


Table 5: Parameters for /3 = 5.40. 


extended the constant rn = mo trajectories down in the direction of the physical 
point. 

While changing the {3 value gives the greatest change to the singlet terms, 
smaller effects occur on the SU{3) flavour symmetric line as we change mo. This 
will help us to locate the initial kq, the starting point for the trajectory in the 
mg-mi plane leading to the physical point. Note that if m is held constant, then a 
further advantage of this condition is that for clover fermions the 0{a) improved 
coupling constant remains unchanged as 

9o = + bgi9o)am), ( 22 ) 

although this is unlikely to lead to any large effect. 

For orientation the SU{3) symmetric point has a pion mass of about ~ 
450 MeV and we reach down to about ~ 260 MeV. 

The specihc components used in the flow discretisation here ar^ 

(f [low], g[auge action], o [bservable]) 

"‘The flow Wilson action here means S'flow[F] = X]FeTr[l — t/P^®-"!]. 
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(W[ilson], S[ymanzik [tree level]], C[lover]) , (23) 


and the Runge-Kutta discretisation is nsed for the flow equation, [2]. 

One can improve the scaling behaviour, which is expected to have O(a^) 
corrections. For example for following |TT] we can write 


F(t) 


1 + + 


= c 


^—^0 imp 


t -t ' 

tOimp — to N- + <^ 2 -^ - -T • • • 

\ -^0 to / 


(24) 


Inverting this and relabelling foimp —>■ Ocont gives 


, _ , ( ^ -^Ocont ® 

to — tocont 1 J- '-'2 

' F/ 


Ocont ^Ocont 


+ . . . 


(25) 


where F^ = F(tx), F^ = tdF{t)/dt\t^, with a; = 0 or Ocont. At tree level for the 
case here {fgo) = (hFS'O), C 2 = —7/72 [18] so we expect the gradient to be +ve. 


4.2 Singlet quantities 

We first investigate the constancy of the singlet quantities X|, as discussed in 
section 121 In Fig. [2] and [3] we plot X| for S = to, wq p and vr against M^/X^ 


0.20 


0.15 


0.10 


0.05 


0.00 


{P,K„) = (5,80,0.122810) 


• S = t, 
► S = N 

■ S = Wo 

◄ S = p 
AS = Jt 


0.0 0.2 0.4 0.6 0.8 1.0 


Figure 2: Top to bottom for S = to (circles), N (right triangles), wq (squares), 

p (left triangles) and tt (up triangles) for (/3,^o) = (5.80,0.122810) with constant fits. 
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0.30 


0.25 


{P,k„)={5.50,0.120900) 


0.20 

r' 0.15 

CO 

X 

0.10 



•s-t. 

► S-N 
■ S = w„ 

•<s.p 

S-n 


A 


0.05 


0.00 


0.0 


0.2 


0.4 0.6 

m:/x: 


0.8 


1.0 


Figure 3: for S = to, N, Wo, p and tt, (top to bottom) for (/?, acq) = 

(5.50,0.120900) together with constant fits. Same notation as for Fig. [2l Opaque 
points are not included in the fits, as they have < 4. 

(which is equivalent to 1/ki). The value at M^/X^ = 1 corresponds to the 517(3) 
symmetric point and the vertical dashed lines correspond to the physical point. 

No structure or trend is seen in the results and they are compatible with 
(X^®'*)^ (and hence Xg^*) a constant down to the vicinity of the physical point. 
The constant fits exclude points with < 4. However these additional points 
all have > 3 and are completely consistent with the fitted points. Thus the 
constancy of X^ as discussed in sections [2] and [3] is supported by the numerical 
results. 

Furthermore any (significant) quadratic term would mean that each quan¬ 
tity S starts at a slightly different point on the SU{3) flavour symmetric line 
(but all with the same gradient) and the trajectories would then all focus at 
the experimental point. We have considered the partially quenched expansion 
(up to cubic terms in the quark mass) and have determined that along the uni¬ 
tary line considered here these higher order terms are negligible - only when 
the quark mass is in the vicinity of the charm quark mass do these non-linear 
terms become appreciable. Thus practically we have a unique starting point on 
the SU{3) flavour symmetric line for the trajectory to the physical point, i.e. 
we have the situation for at least S = to, N, wo, p and vr of the left panel of 
Fig. [U To check this an alternative description is provided by a re-arrangement 
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(26) 


of XUXI = {2Ml + Ml)/XI to give 

2MI-MI XI ^Ml 

Ag Ag Ag 

for S = N, p,to,Wo. So plotting {2M]^ — MD/X^ against Ml/X^ with constant 
gradient —2 should describe the data for all S. Hence in Fig. [T] left panel the 
gradient is —2 while for the right panel, we would have an initial gradient of —2 
and then some curvature, but all meeting at the physical point. 

In Figs. [Hand [5] we present hts of eq. fl26|l to numerical results of (2M]^ — 



M//X/ 


Figure 4: {2Ml^ — Ml)/Xg against Ml/Xg, together with the fit from eq. ([2U|1 for 
(/3,ko) = (5.80,0.122810) and S = N (right triangles), p (left triangles), to (circles) 
and Wo (squares). The stars correspond to the experimental values for S = p, N upper 
and lower respectively. 

M^)/X| versus Ml/X^ with S = N, p, to, Wo for the same data sets as in Fig. [2] 
and |3l Again straight lines (from the £t function) describe the data very well. 
However as can be seen the lines for the S = N and p cases do not quite go 
through their physical points (denoted by stars). This is because the kq used 
while close is not quite the value required for the correct path. We shall in 
future denote this point by Kq. For example, it can be seen that the [3 = 5.80, 
Kq = 0.12281 lines are closer to the physical point, i.e. Kq = 0.122810 is closer to 
Hq than for the (3 = 5.50, kq = 0.120900 data. Again the picture is best described 
by the left panel of Fig. [H So at least for all these quantities we only have to 
slightly tune to find the appropriate Kq giving the beginning of the path from the 
SU{3) flavour symmetric line to the physical point. 
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m//x/ 


Figure 5: (2M|^ — M‘^)lX‘g against M^/Xg, together with the fit from eq. ([26]l for 
{/3,ko) = (5.50,0.120900) and S = N, p, to, wq. The notation is as for Fig. [H Opaque 
points are not included in the fit as they have < 4. 

5 Scale setting 

Using these results we now take Xs = const, to define the scale, i.e. we set 

= const. = asXg"^ . (27) 

We take the experimental hadron mass results as given in Table [H As Xtg, 
are secondary quantities, i.e. X“p are not experimentally known, they have 
to be determined. 

If we now normalise 

®s(^o) / .^exp\2 ; 

) 

this provides an estimate for the lattice spacing using the singlet quantity S, 
which is also a function of the point on the SU{3) flavour symmetric line, i.e. kq 
(we have indicated this by writing 05 (^ 0 ) for the lattice spacing). 

We now vary kq, searching for the location where the various 05 (^ 0 ) cross, 
providing a value for the common lattice spacing a (and Kq). While ideally we 
would wish the crossing of all the lines to occur at a single point leading to a 
common lattice spacing, this, of course, does not quite happen. So we consider 
pairs of singlet quantities and determine the crossing points, together with the 
associated (bootstrap) error. In particular we apply this to the pairs 

(7r,X), (7r,p) . (29) 
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We now use these crossings to adjust Xt^ and so they also go through 
these points. This determines For example we have 


[w. 


exp\2 _ 
0 ) = 




0.004 


II 

^ CO 

03 


0.003 


0.002 


8.138 8.140 


(30) 


For example in Figs. Eland [7] we plot a|(/Do) from eq. fl25]) (in fm^) against 

0.005 


A 

|3=5.80 (ji.N) 1 

\ ... - - 

._: 

A 

► S = N 

S = Jt 
• s = to 
■ S = Wo 


8.142 8.144 

1/k„ 


8.146 8.148 


Figure 6: a| against I/kq for 5 = tt (up triangles), N (right triangles) and to (circles), 
wq (squares) together with linear fits for /3 = 5.80. 

kq for {3 = 5.80 and 5.50 for the singlet quantities S = n and N together with 
S = to and wq- Where we have three kq values a linear £t in I/kq is made, while if 
there are four ko values available then a quadratic £t is made. (However it made 
very little difference to the later results whether the results from the linear or 
quadratic fit is used, as mainly interpolations between the data is sufficient.) 
Also plotted is S' = to and Wq, again together with appropriate fits. The lattice 
values have been adjusted with a common factor so these singlet quantities also 
cross at the same value as (tt, N), which is equivalent to a determination of 
and Wo'^'^ as indicated in eq. fl5U]l . This procedure is then repeated for the pair 

(tp)- 

For completeness we also take a weighted average of both the (tt, N) and (tt, p) 
crossings to determine the best (l/«)g,a^). These values are given in Table El 
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0.007 


0.006 


'i” 0.005 


0.004 



8.264 8.266 


8.272 


Figure 7: Og against I/kq for S = tt, N and to, wq together with quadratic fits for 
P = 5.50. Notation as for Fig. [6l 



1 /k*, 

[fm"*] 


a [fm] 

5.80 

5.65 

5.50 

5.40 

8.14197(12) 

8.19602(15) 

8.26844(13) 

8.33823(25) 

0.00346(04) 

0.00468(06) 

0.00547(06) 

0.00669(16) 

0.122820(2) 

0 .122010(2) 

0.120942(2) 

0.119930(4) 

0.0588(03) 

0.0684(04) 

0.0740(04) 

0.0818(09) 


Table 6: Determined values of I/kq a? [fm^]. For completeness in the third and 
fourth columns we also give Kq and a [fm] directly. 

6 Continuum results 

We are now in a position to perform the last, continuum, extrapolation. In 
Figs. [HI E] we show these extrapolations from the pairs (vr,iV) and (7r,p). As 
anticipated the gradients in a?, while small, are positive (c.f. eq. fl25|) i with the 
(tt, p) results being slightly larger than the (vr, N) results. 

Finally a weighted average from these continuum results (i.e. for = 0) gives 
our hnal results 

= 0.1511(22)(06)(05)(03)fm, = 0-1808(23)(05)(06)(04) fm. (31) 

The hrst error is statistical, while the second (hnite volume), the third {SU{3) 
flavour breaking expansion) and fourth (scale) are systematic errors as discussed 
in Appendix [Bl 
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Figure 8: y/to and ioq (in fm) against (in fm^) from the {tt,N) crossing together 
with a linear fit. 
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Figure 9: ^/tQ and wq (in fm) against (in fm^) from the (vr,/?) crossing together 
with a linear fit. 
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7 Conclusions 


In this article we have described a method for determining the trajectory to 
approach the physical point, and demonstrated (theoretically and numerically) 
that singlet quantities remain constant as we approach this point. This enables 
us by considering pairs of singlet quantities to determine ‘best’ lattice spacings 
and starting values for the path. By matching these results to the flow variables 
to and tco this enables a determination of their physical values, see eq. (1^ . 

In Fig. [10] we compare these results with other determinations for Uf = 2 + 1 


0.14 


BMW 12 


HotQCD 14 


RBC-UKQCD 14 


QCDSF-UKQCD 15 


0.15 

[fm] 


0.16 0.17 


0.18 
w„ [fm] 


0.19 


Figure 10: left panel and right panel in fm for BMW 12 [3], HotQCD 

14 [9j, RBC-UKQCD 14 [TO], together with the present results. 

flavours, namely BMW 12 [3], HotQCD 14 [9] and RBC-UKQCD 14 (The 
given errors are taken in quadrature.) Reasonable consistency is found between 
the different determinations. 

In conclusion we have determined in this article the flow scales for to and Wq. 
These are ‘secondary’ scales and while having the advantage of being cheap and 
accurate to determine from lattice simulations are not directly experimentally 
accessible and thus have to be matched to physical quantities. 


Acknowledgements 

The numerical conhguration generation (using the BQCD lattice QCD program 
[23]) and data analysis (using the Chroma software library [2l]) was carried out 


17 









on the IBM BlueGene/Q using DIRAC 2 resources (EPCC, Edinburgh, UK), the 
BIueGene/P and Q at NIC (Jiilich, Germany), the Lomonosov at MSU (Moscow, 
Russia) and the SGI ICE 8200 and Cray XC30 at HLRN (The North-German Su¬ 
percomputer Alliance) and on the NCI National Facility in Canberra, Australia 
(supported by the Australian Commonwealth Government). HP was supported 
by DEG Grant No. SCHI 422/10-1. PEER was supported in part by the STFC 
under contract ST/G00062X/1 and JMZ was supported by the Australian Re¬ 
search Council Grant No. FT100100005 and DP140103067. We thank all funding 
agencies. 

Appendix 

A Singlet chiral perturbation theory: Wilson 
flow 

We want to check that the chiral perturbation theory results for the Wilson flow, 
as given in 121 ] are consistent with the SU{3) flavour symmetry expansion [I6] . 

A.l Pseudoscalar Meson masses 

First we need to set out some notation. The quark masses for the 2-1-1 case are 
best denoted by dehned through 


Xi = Bomi 
Xs = BoTHs . 


(32) 


To simplify expressions, it is useful to dehne some additional y variables: 


X = |(2Xi + Xs) 

Xn = Xl 


Xk = ^{Xs + Xi) 

Xv = \{‘^Xs + Xi), 


(33) 


and a logarithmic function 



P e n, K^r]. 


( 34 ) 


In this notation the NLO pseudoscalar meson masses are [22] 


Ml 
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where 


Mk = XK\l + qlX + q2XK + ‘^^^■n^ 

= Xr?|l + giX + ?2Xr? + 2/iit- 

+ + ‘^l^K + ^/ir,| + qsiXs - Xif 


qi 

q2 

q3 


—^{2Le — L 4 ) 
JO 
16 

71 ( 2 ^ 8 - 1 ^ 5 ) 

Jo 

772(8-^7 + Ls) 

yjo 


( 35 ) 


(36) 


A.2 Wilson Flow scale, to 

In [21], eq.(4.10), Bar and Golterman give the form expected for the quantity to 
at NNLO in chiral perturbation theory. 


tn — t, 


0,ch 


1 + 


ki 


(4vr/o) 


-X2MI + MI 


+ ^(3^2 — + 4:k2M]^fiK + ~ ^MK)jl'rj + ^2^^/7r? 


*=* (2MJ + m;)= + j-Xr^iMl - Ml)^ 


(4vr/o) 


(4vr/o)^ 


(37) 


The free parameters in this expression are ki,k 2 ,ki,k^. Most terms in the ex¬ 
pression are obviously symmetric, but the ki term has been written in a way that 
obscures its symmetry. 

Using eqs. fl35ll to translate eq. flTTl) into x variable^ 


tn — t, 


O.ch 


3ki _ ^2 /o \ 

1 H ~ TTTr X { A — £ ^9 v^XTr/^TT “1“ ^XKf^K H“ Xr]f^7]) 


{Anfor " ( 477 / 0 )^ 


9 /C 4 


(4vr/o) 




ki 


4(47r/o)^ 


(Xs - XiY 


with 


k'^ = k4 + liAnfofkiiqi + ga) 

^5 = ^5 + |(4vr/o)^/cig2 • 

®In this Appendix we work to order Xq, dropping terms of order Xq- 


(38) 


( 39 ) 
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The expression in eq. fl38D is simpler and more explicitly symmetric than eq. fl37|) . 

As usual, we get a further simplihcation if we restrict ourselves to the line of 
constant y, 


with 


and 


tn — T 


1 + 


3x1 In ^ + 4:Xk In ^ In ^ 


(47r/o)^ V ^ T 

u> 

^ iXs - Xif 


X 


X 


4(47r/o)^ 


k'l = kl + ^^2 In ^ 


9 


A2 


T — to, 


ch 


1 + 


3^1 _ , 8k2 _2, X , 

■X + Tir^XlX Iny^T 


9k A 


(4vr/o)2 (47r/o)^ 


A2 (47r/o)4 




(40) 


(41) 


(42) 


being the value of to on the symmetric line. We can Taylor expand eq. fl40|) about 
the symmetric point, the result is 


to 


T 


1 + 


(4vr/o) 


{lk2 + \k'i){xs-Xir + --- 


4^6 


(43) 


As expected, there is no linear term, and the first term we see is quadratic in the 
SU{3) breaking. 


B Systematic errors 

We follow here the more general discussion given in Appendix A of |19] . 

B.l Finite lattice volume 

Clearly the argument given in section [2] that Xs is flat along the SU{3) flavour 
symmetric point holds for any volume. As discussed in [16] for an estimate of 
hnite volume effects, a suitable expression is given by 

XHL) = A'l (l + Csi[/i(M,) + 2fdMK)]) . (44) 

Lowest order xPT; |25l |26] indicates that reasonable functional forms for fiiM) 
are 

^-ML 

/l(M) = {aMf , meson, 

e-ML 

hiM) = {aMf , baryon. 
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(45) 
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Figure 11: versus (/^(MTr) + 2/2,(Mx))/3 for (/?, kq) = (5.50,0.120900, with 

N (squares), p (diamonds) and tt (upper triangles). The left-most clusters of points 
are from the 32^ x 64 and 48^ x 96 lattices, while the right cluster is from 24^ x 48 
lattices. The dashed lines are linear fits. 


In Fig. [H] we plot {fiiM.,,) + 2/i(Mi^))/3 against (X^*)^ for S = N, p and vr on 
48^x96, 32^x64 and additionally 24^x48 lattices for /? = 5.50 and kq = 0.120900. 
The fits are linear, with reasonable agreement to the data. Little hnite size effect 
is seen between the larger lattice volume results used in the previous analysis and 
the extrapolated value here. For the case considered here for {X.„, X^, Xp) the 
changes are about (1.3, 2.9,1.5)%. Taking this as an increase in errors for all the 
data sets and performing the same analysis gives the change in the central values 
and errors, taken as the systematic error as given in eq. fl3Tl) . 

B.2 SU{3) flavour breaking expansion 

We hrst note that in Figs. [2, [3] from the SU{3) flavour symmetric line down to 
the physical point lies in the range \6mi\ < 0.01 (and \6ms\ < 0.02), e.g. [20] 
and that mass ‘fan plots’ (e.g. Fig. 5 of [20|) show little curvature. This is in 
agreement with the SU{3) flavour breaking expansion, eq. ffT^ or ffTOD . The next 
order in the expansion is multiplied by a further Srriq. So we expect that every 
increase in the order leads to a decrease by an order of magnitude or more (often 
by a factor ~ 20) in the series. So we believe that convergence is very good for 
hyperons. (Such an expansion is good compared to most approaches available 
to QCD.) Nevetheless we have, however, made tests with a linear or quadratic 
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fit for example for the nucleon in Fig. |3]and followed this through the analysis. 
The hnal change in central value for y/to and wo was not large, we include it as 
a (second) systematic error. 

B.3 Physical scale 

As mentioned in footnote |2] physical values of hadron masses have a small elec- 
tromagetic component. Although we disregard this in our analysis, we make a 
small allowance here, and take to also have a similar error as i.e. 

a systematic error of ~ 0.2% due to electromagnetic effects. 
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